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SUMMARY 


Ice accretion calculations were made for a modem commercial transport using the NASA, 
Lewis LEWICE3D ice accretion code. The ice accretion calculations were made for the wing and 
horizontal tail using both isolated flow models and flow models incorporating the entire airplane. 
The isolated flow model calculations were made to assess the validity of using these simplified 
models in lieu of the entire model in the ice accretion analysis of full aircraft. Ice shapes typifying 
a rime and a mixed ice shape were generated for a 30 minute hold condition. In general, the calcu- 
lated ice shapes looked reasonable and appeared representative of a rime and a mixed ice condi- 
tions. The isolated flow model simplification was good for the main wing except at the root where 
it overpredicted the amount of accreted ice relative to the full aircraft flow model. For the horizon- 
tal tail the size and amount of predicted ice compared well for the two flow models, but the posi- 
tion of the accretions were more towards the upper surface for the aircraft flow model relative to 
the isolated flow model. This was attributed to downwash from the main wing which resulted in a 
lower effective angle-of- attack for the aircraft horizontal tail relative to the isolated horizontal 
tail. 


NOMENCLATURE 


Airplane angle-of-attack, degrees 
Span of wing, m 
Chord length, m 
Pressure coefficient 
Droplet diameter, pm 
Roughness factor, m 
Liquid Water Content, g/m 3 
Median Volume Diameter, pm 
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P Free stream pressure, Pa 

s Surface distance, cm 

t Ice accretion time, seconds 

T Free stream temperature, K 

V Free stream velocity, m/s 

y Spanwise distance from root, m 

(3 Collection efficiency 


I. INTRODUCTION 

The problem of aircraft icing has plagued man since the beginning of flight. Ice, in general 
causes increased drag, decreased lift, decreased propulsive thrust, increased weight, reduced sta- 
bility, etc., all of which can make an aircraft less controllable and more prone to fall out of the sky. 
Simultaneously with the problem of icing came the desire to rid ourselves of this problem, which 
led to the task of trying to understanding the ice accretion phenomena. To better protect aircraft 
from the icing, the icing phenomena needed to be understood. Early attempts at understanding the 
cause and effects of icing included various flight tests and some open air icing tests atop Mount 
Washington (Ref. 1). Some early analytical work was done by Langmuir and Blodgett to calculate 
the collection efficiencies of cylinders using a mechanical integrator with numerous simplifica- 
tions (Ref. 2). Much was done to quantify the effect of icing, but limited insight was gained into 
the underlying physics of icing due to deficiencies in the analytical and experimental tools avail- 
able to the early researchers. From this early work, conservative engineering tools were devel- 
oped which allowed design of ice protection systems and allowed aircraft to fly safely. 

With the advent of high speed computers in the past sevaral decades and the desire to opti- 
mize system design, the task of analytically predicting the icing phenomena has been undertaken. 
Many computer codes have been written in the past decade which calculate ice accretions for 2D 
airfoils (Ref. 3-6). Calculations of ice accretion for some simple 3D geometries have been made 
in the past several years (Ref. 7-9). These cases have been limited to isolated portions of an air- 
craft and mainly to the wing and tail sections. The work presented here is the first work which 
attempts to analyze the entire aircraft and the first work which attempts to assess the validity of 
the isolated model simplification. Ice accretion predictions on the wing and horizontal tail were 
made using the aerodynamic model of the entire aircraft and using isolated aerodynamic models. 
A typical commercial transport configuration was used for the calculations. 

n. ANALYTICAL METHOD 

The NASA, Lewis ice accretion code (LEWICE3D) was used to do the ice accretion anal- 
ysis (Ref. 10). This computer program has been used in previous calculations of isolated, finite 
wings (Ref. 7-9). This work represents the first application of the code to the full aircraft. 

The LEWICE3D code incorporates flow, trajectory, heat transfer and geometry modifica- 
tion analysis into a single computer program with a multitude of capabilities. Flow can be calcu- 
lated for arbitrary 3D lifting and nonlifting bodies with external flow. Calculations of arbitrary 
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streamlines and trajectories are possible. The code has the capability to calculate tangent trajecto- 
ries and impingement efficiencies for single droplets or droplet distributions. Ice accretions can be 
calculated at arbitrary regions of interest for lifting surfaces in a multistep fashion in either a sur- 
face normal or tangent trajectory direction. Additionally, geometry for lifting surfaces (i.e. airfoil 
surfaces) can be modified automatically at each time step to reflect the ice shape. 

The methodology used in the LEWICE3D analysis can be broken into seven basic steps 
for each section of interest at each time step. In the first step the flow field is generated. Secondly, 
surface streamlines are calculated. Thirdly, tangent trajectories are calculated at the region of 
interest. An array of particles is released between the tangent trajectories in the fourth step. These 
impacting particles are used to calculate collection efficiency as a function of surface distance. 
The fifth step involves interpolating or extrapolating the collection efficiencies onto the stream- 
lines. In the sixth step the ice accretion for the streamline is calculated. During the seventh and 
final step the geometry is modified. 

There are four basic program elements contained in the LEWICE3D jobstream; a flow 
analysis, a trajectory analysis, a streamline analysis and an ice accretion analysis. The flow code is 
basically the 1st order, 3D panel code of Hess-Smith (Ref. 1 1-12) for arbitrary lifting and nonlift- 
ing bodies with external flow with modifications by Hillyer Norment (Ref. 14) to allow leaking 
panels (used to simulate mass flows through orifices). It incorporates variable dimensioning to 
allow easy adoption to different sized computers and problems. The trajectory analysis is basi- 
cally that of Hillyer Norment with modifications by Bidwell and Mohler (Ref. 10) to the velocity 
calculation (vectorized) and by Bidwell to the tangent trajectory search scheme (replaced with 
more efficient scheme). At the heart of the trajectory analysis is the variable step predictor-correc- 
tor integration scheme by Krogh (Ref. 15). The surface streamline analysis uses a variable step 
size fourth-order Runge-Kutta integration scheme developed by Bidwell (Ref. 10). The ice accre- 
tion model is basically that of the LEWICE2D code applied along surface streamlines (Ref. 4, 

16). 

m. CONFIGURATION 

The conditions and geometry for the analysis were chosen to typify a commercial trans- 
port in a hold condition. The configuration used in the analysis, which was similar to a Boeing 
737 airplane, was chosen to be consistent with those of a wind tunnel model to be used at Langley 
Research Center for iced stability and control measurements. The coordinates for this model were 
obtained from Simha Dodbele at Langley Research Center. A hold condition was chosen because 
it was a crucial design point for the ice protection systems on a commercial transport. 

The aircraft model used was a simplified version of the Boeing 737 airplane. The engine 
nacelles and engine pylons were not included in the flow model. The wings were modeled using a 
single element model with retracted flaps and slats. The deletion of these features from the aircraft 
should not be deemed too severe considering the position of the engines and the position of the 
flaps and slats in a typical hold condition. 

Two flight conditions were chosen to explore limiting points of interest in the icing enve- 
lope for a hold condition. These points yield a typical rime and mixed icing condition. For both 
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flight conditions the airplane angle-of-attack, yaw angle, altitude, and airspeed were chosen to be 
0 degrees, 0 degrees, 1000 m, and 135 m/s respectively. The hold or icing time for both conditions 
was chosen to be 30 minutes. In addition, for both flight conditions a Median Volume Diameter 
(MVD) of 20 microns was chosen as indicative of a typical drop size. A temperature of 243.1 K 
and an LWC of .20 g/m3 were chosen for the rime case. A temperature of 263.7 K and an LWC of 
.695 g/m3 were chosen for the mixed case. 

The computer program parameters were chosen from experience, correlations and a desire 
to limit the computational resources required. A single drop size and icing time step were chosen 
for the calculations The calculations were made at 6 spanwise stations on the wing and at 3 span- 
wise stations on the horizontal tail. In the wing analysis 3 spanwise stations were distributed on 
the outer portion of the wing and 3 were distributed on the shoulder or inner portion of the wing. 
The spanwise stations were chosen to be at the 10%, 50% and 90% stations on the outer wing, the 
inner wing, and the horizontal tail. The 10% and 90% span stations were chosen as conservative 
spanwise limits for the LEWICE3D methodology. Outside of these limits spanwise pressure gra- 
dients become large and the strip theory assumptions built into the LEWICE3D ice accretion code 
can be violated. The full aircraft calculations for the wing and horizontal tail were done separately 
using different aerodynamic models. The aerodynamic models used for each of the cases was 
refined in the area where the ice accretions were to be calculated. This allowed for fewer panels 
and smaller computational times. 

Additional calculations were made using isolated models of the wing and tail. These were 
done to understand the correlation between the full airplane result and the isolated airfoil result. 
The isolated airfoil simplification is commonly used in the icing analysis of aircraft to reduce cost 
and complexity. A typical isolated model may have 60% fewer panels than its full aircraft coun- 
terpart, resulting in a 60% savings in computer time. The calculations were made on the NASA 
Lewis Research Center CRAY XMP and on an IRIS Model 4D/440/VGX. The calculations for 
the entire aircraft, which included 6 stations on the wing and 3 stations on the horizontal tail, 
required approximately 10 hours of CPU time on the Cray XMP and approximately 100 hours of 
CPU time on the IRIS workstation. The isolated calculations, which were made at the same loca- 
tions, required approximately 4 hours and 40 hours of CPU time on the CRAY XMP and IRIS 
workstation respectively. 

IV. ANALYSIS 

Two icing conditions and two aerodynamic flow models were used in the analysis of the 
wing and horizontal tail. The two icing conditions were chosen to explore a rime and a mixed ice 
hold condition for the wing and tail. Two flow models were used in the analysis of the wing and 
tail elements; an isolated model of the element (“isolated” model), and a model with the element 
and the remainder of the aircraft modeled (the “aircraft” model). The isolated flow models were 
used to assess the validity of using the isolated flow model instead of the actual flow model in the 
icing analysis. 

The wing analysis panel models are shown in Figures 1 and 2. The isolated model con- 
tained 4320 panels while that of the aircraft model contained 7292 panels. Figure 3 depicts the air- 
foil section at the five spanwise stations where the analysis was done. In general the wing was 
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complex, having taper, twist and airfoil section variation. Figure 4 shows the pressure distribu- 
tions at each of the spanwise stations for the models. For both models the section lift coefficient 
increases as the root is approached due to the twist in the wing. The agreement in the pressure dis- 
tribution at the tip is almost perfect for the two models, but falls off as the root is approached. The 
section lift coefficient for the isolated model is greater than for the aircraft model due to the effect 
of the fuselage. 

The collection efficiency for the two aerodynamic wing models is shown if Figure 5. In 
general, for both models, the shape of the collection efficiency curves were similar except for the 
most inboard section. The curves show that the majority of water is collecting on the upper sur- 
face (i.e. a negative surface distance from the highlight) indicative of a negative angle-of-attack, 
although the wing was at a slight positive angle-of-attack. This was due to the low angle-of-attack 
and camber of the wing. The shapes of the collection efficiency curves are similar at the outer four 
stations because the airfoil sections are similar. The airfoil sections at the two inboard sections are 
much blunter hence the collection efficiency curves are more rounded. In general, for the three 
outer stations, the maximum collection efficiency decreased while the upper and lower impinge- 
ment limits moved more towards the underside or pressure side of the wing as the root was 
approached. This is due to the increased inertia parameter (resulting from the increased chord) 
and effective angle-of-attack as the root is approached. At the fourth station, which has the same 
section as the outer three stations, the maximum collection efficiency decreased relative to the 
outer three station but the extent of impingement increased slightly. This can be attributed to the 
increase in the sweep angle of the leading edge at the inboard portion of the wing. Increased 
sweep angle causes a reduction in maximum collection efficiency and an increase in extent of 
impingement. At the fifth and sixth stations the maximum collection decreased and the extent of 
impingement increased over the outer stations. This decrease in maximum collection efficiency 
can be attributed to the increase in chord, while the increased extent of impingement can be attrib- 
uted to the increased sweep angle and the blunter section shape. 

As was the case for the pressure distribution, the agreement in collection efficiency 
between the two models is almost perfect at the tip, but worsens as the root is approached (Figure 
6). Only at the inner two stations is there in appreciable difference in the collection efficiency 
curves. At the fifth station for both models, the maximum collection efficiency and the lower 
impingement limits are about the same, while the upper impingement limit of the airplane wing is 
increased slightly over the lower impingement limit of the isolated wing. This is due to a differ- 
ence in effective angle-of-attack between the two models. The airplane wing at this station has a 
slightly smaller effective angle-of-attack than the isolated wing hence the droplets impinge further 
back on the upper surface. At the sixth and most inboard station the maximum collection effi- 
ciency and upper impingement for both models is about the same while the lower impingement 
limit is further back for the isolated model. This can be explained once again by the fact that the 
aircraft wing is at a lower effective angle-of-attack than the isolated wing. The difference in the 
way that the lower angle-of-attack is revealed in the collection efficiency curves between the fifth 
and sixth section, (i. e. one shows a much different upper impingement limit and one shows a 
much different lower impingement limit) is due to the radically different section shapes. The most 
inboard section is much more rounded on its underside hence a small decrease in angle-of-attack 
at a relatively low angle-of-attack will produce a larger change in the lower impingement limit. At 
the fifth station the section is less rounded on its underside, hence a small decrease in effective 
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angle-of-attack at a relatively low angle-of-attack will show up more upon the upper impingement 
limit. 


The ice accretion shapes for both the rime and mixed conditions and for both flow models 
are shown in Figures 7-10. For both the rime and the mixed conditions a roughness factor of .0013 
m was used in the ice accretion calculations. This value had been used in previous calculations 
with good agrement (Ref. 8). Doubling or halving the roughness factor did not significantly affect 
the resulting ice shape, and in light of the absence of experimental data, and the reasonable 
appearance of the resulting ice shapes it was deemed a good value. In general, the rime and mixed 
shapes looked reasonable considering the collection efficiencies and the pressure distributions 
observed. The ice shapes followed the trends observed for the collection efficiency. That is, as the 
root was approached, the ice shape height decreased and its extent increased or decreased accord- 
ing to the local extent of impingement. The comparison between the isolated and aircraft wing ice 
shapes was excellent at all but the root section. This agreement was expected considering the 
good agreement between collection efficiency and pressure distribution for the isolated and air- 
craft wing models. At the root the differences in the pressure distribution and collection efficiency 
produce some differences in the resulting ice shape between the isolated and aircraft wing, 
although not much. 

The flow models for the isolated and aircraft horizontal tail are shown in Figures 1 1 and 
12. The isolated panel model contained 1800 panels while that of the aircraft model contained 
5532 panels. The airfoil sections for the three spanwise sections analyzed are shown in Figure 1 3. 
The horizontal tail was of a simple tapered design. The pressure distributions for both models are 
shown in Figure 14. The pressure distributions reveal a negative lift coefficient. This is due to the 
relatively low aircraft angle-of-attack, the inverse camber of the tail section (i.e. the tail camber 
was in the opposite direction of the wing camber), and the fact that the horizontal tail is mounted 
at a slightly negative angle-of-attack relative to the aircraft angle-of-attack (about 1 degree). Also, 
the pressure coefficients display a slight fall off in section lift as the root is approached, a result 
indicative of a simple tapered wing with no twist. From the pressure distributions it can be seen 
that the effective angle-of-attack for the isolated model is less than that for the actual model. The 
difference in effective angle-of-attack increases as the root is approached. This difference is 
caused by downwash from the main wing onto the horizontal tail. 

The collection efficiencies for the two tail models are shown in Figure 15. In general, the 
maximum collection efficiency decreased, the extent of impingement increased, and the region of 
impingement moved more towards the upper surface of the wing as the root was approached. The 
decrease in the maximum collection efficiency and increase in the extent of impingement is due to 
the increase in chord length as the root is approached. The collection efficiency for the horizontal 
tail, which is relatively blunt and at a low angle-of-attack, behaves much like a swept, tapered cyl- 
inder with increasing diameter towards the root. The maximum collection efficiency for a tapered 
swept cylinder decreases and the extent of impingement increases as the root is approached. The 
migration of the impingement region towards the upper surface as the root was approached is due 
to the spanwise decrease in section lift or effective angle-of-attack. 

As was the case for the wing, the collection efficiency differences between the isolated 
and aircraft flow model for the horizontal tail were small at the tip but increased as the root was 
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approached (Fig. 16). As the root is approached the difference in maximum collection efficiency 
between the two models increased, with the actual model producing lower values. Also as the root 
was approached the region of impingement for the aircraft model moved increasingly towards the 
upper surface of the airfoil relative to the isolated model. These differences can once again be 
attributed to the differences in effective angle-of-attack. The aircraft model saw a smaller effec- 
tive angle-of-attack than the isolated model and this difference in effective angle-of-attack 
increased as the root was approached. Hence the region of impingement of the aircraft model rel- 
ative to the isolated model was more towards the upper surface and it moved increasingly towards 
the upper surface as the root was approached 

The ice accretion shapes for both the rime and mixed conditions for both horizontal tail 
flow models are shown in Figures 17-20. As for the wing ice accretions, a roughness factor of 
.0013 m was used to produce the ice shapes. The ice accretions, appeared reasonable and repre- 
sentative of rime and mixed conditions and followed the trends shown by the collection efficien- 
cies. That is, the ice shape size decreased and extent of impingement increased as the root was 
approached. The agreement in size and shape between the ice accretions for the isolated and air- 
craft models was good. The ice accretions for the actual model were more towards the upper sur- 
face relative to the isolated flow model due to the lower effective angle-of-attack of the aircraft 
model relative to the isolated model. 

V. CONCLUSION 

In general, the calculated results were encouraging. The calculated flow looked reasonable 
considering the angles-of-attack and the potential flow assumptions used. The calculated collec- 
tion efficiencies and ice shapes were consistent with previous work and intuition. Comparisons 
between the isolated flow model and the aircraft model were good for all but the most inboard sta- 
tions for the wing and horizontal tail. 

The flow solutions produced by the panel code appeared reasonable and followed tradi- 
tional trends attributed to finite swept, tapered, twisted wing. The majority of the horizontal tail 
and wing were of the same cross section, hence variation in the pressure distribution could be 
attributed to effective angle-of-attack differences caused by geometric twist, taper and down wash. 
In all cases the size and type of variation observed could be attributed to one of these effects. 

Several trends with regard to collection efficiency were observed in the study. Because the 
airfoil sections and sweep angle for a majority of the wing and tail were similar the trends in col- 
lection efficiency could be mainly attributed to chord length, and local angle-of-attack. The maxi- 
mum collection efficiencies decreased as the root was approached for both the wing and tail 
models due to the increase in chord length. The limits of impingement for the swept, tapered wing 
could be correlated to the local effective angle-of-attack. The effective angle-of-attack of the wing 
increased as the root was approached due to the geometric twist, hence the limits of impingement 
moved more towards the pressure side of the wing. The extent of impingement on the horizontal 
tail, which was blunt and at a low angle-of-attack, increased as the root was approached due to the 
increase in chord. 

In general the calculated ice shapes were reasonable and representative of the rime and 
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mixed conditions from which they were derived. The ice shape height or maximum thickness 
decreased as the root was approached for both conditions following the maximum collection effi- 
ciency trend. Extent of ice shape also followed the extent of impingement trend. 

The pressure distribution, collection efficiency and ice accretion comparisons between the 
isolated and actual flow models were good at all but the most inboard sections where the isolated 
flow model yielded more conservative results. The differences in the pressure distribution, collec- 
tion efficiency, and ice accretions between the isolated and aircraft sections can be attributed to 
affective angle-of-attack differences between the models. In the case of the wing, the fuselage 
produced a slightly lower effective angle-of-attack relative to the isolated wing. For the horizontal 
tail the downwash from the wing caused smaller effective angle-of-attacks relative to the isolated 
tail. The quality of the isolated calculations were good relative to the aircraft calculations for the 
work presented, but caution should be used in extrapolating these results to higher angles-of- 
attack. At the higher angles-of-attack the differences in the effective angle-of-attack for both mod- 
els will increase resulting in poorer comparisons. Methods to correct the effective angle-of-attack 
in the isolated cases need to be developed to make the isolated simplification more viable. 
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PRESSURE COEFFICIENT (CP) PRESSURE COEFFICIENT (CP) 



(a) AIRCRAFT WING. 



(b) ISOLATED WING. 


FIGURE 4. - PRESSURE DISTRIBUTION AT SEVERAL SPANWISE LOCATIONS ON WING FOR AAOA 
0 DEGREES. 
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(a) AIRCRAFT WING. 



SURFACE DISTANCE FROM HIGHLIGHT (CM) 


(b) ISOLATED WING. 

FIGURE 5. - COLLECTION EFFICIENCY AS A FUNCTION OF SURFACE DISTANCE AT SEVERAL 
SPANWISE LOCATIONS ON WING. FLIGHT CONDITIONS; AIRSPEED. 135 M/S, AAOA; 0 DEGREES, 
DROP SIZE; 20 urn, STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa. 
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FIGURE 6. - COMPARISON OF COLLECTION EFFICIENCY FOR THE ISOLATED AND AIRPLANE 
WING AS A FUNCTION OF SURFACE DISTANCE AT SEVERAL SPANWISE LOCATIONS ON WING. 
FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 20 jim, STATIC TEM- 
PERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa. 


SPANWISE STATION (y/b) 



FIGURE 7. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON AIRCRAFT 
WING. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 20(j.m, STATIC 
TEMPERATURE; -29.9 C, STATIC PRESSURE; 89867 Pa. ICING TIME; 30 MINUTES. LWC; .2 G/M 3 . 
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SPANWISE STATION (y/b): 



FIGURE 8. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON ISOLATED 
WING. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, A AO A; 0 DEGREES, DROP SIZE; 20 4m, 
STATIC TEMPERATURE; -29.9 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. LWC; 
.2 G/M 3 . 


SPANWISE STATION (y/b): 



FIGURE 9. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON AIRCRAFT 
WING. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, A AO A; 0 DEGREES, DROP SIZE; 20 4m, 
STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. LWC; 
.695 G/M 3 . 
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FIGURE 10. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON ISOLATED 
WING. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 20 urn, 
STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. LWC; 
.695 G/M 3 . 


FIGURE 11. - PANEL REPRESENTATION OF AIRCRAFT USED IN HORIZONTAL TAIL CALCULA- 
TIONS. 
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PRESSURE COEFFICIENT (CP) PRESSURE COEFFICIENT (CP) 



(a) AIRCRAFT HORIZONTAL TAIL. 



(b) ISOLATED HORIZONTAL TAIL. 


FIGURE 14. - PRESSURE DISTRIBUTION AT SEVERAL SPANWISE LOCATIONS ON HORIZONTAL 
TAIL FOR AAOA OF 0 DEGREES. 



(a) AIRCRAFT HORIZONTAL TAIL. 



(b) ISOLATED HORIZONTAL TAIL. 

FIGURE 15. - COLLECTION EFFICIENCY AS A FUNCTION OF SURFACE DISTANCE AT SEVERAL 
SPANWISE LOCATIONS ON HORIZONTAL TAIL. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 
0 DEGREES, DROP SIZE; 20 jim, STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa. 
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FIGURE 16. - COMPARISON OF COLLECTION EFFICIENCY FOR THE ISOLATED AND AIRPLANE 
HORIZONTAL TAIL AS A FUNCTION OF SURFACE DISTANCE AT SEVERAL SPANWISE LOCA- 
TIONS ON WING. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 20 
|im, STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa. 



FIGURE 17. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON AIRCRAFT 
HORIZONTAL TAIL. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 
20 urn, STATIC TEMPERATURE; -29.9 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. 
LWC; .2 G/M 3 . 
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FIGURE 18. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON ISOLATED 
HORIZONTAL TAIL. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 
20 urn, STATIC TEMPERATURE; -29.9 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. 
LWC; .2 G/M 3 . 


SPANWISE STATION (y/b): 
— .10 



FIGURE 19. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON AIRCRAFT 
HORIZONTAL TAIL. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, AAOA; 0 DEGREES, DROP SIZE; 
20 \im, STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. 
LWC; .695 G/M 3 . 
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SPANWISE STATION f y/b) : 



FIGURE 20. - CALCULATED ICE ACCRETIONS AT SEVERAL SPANWISE LOCATIONS ON ISOLATED 
HORIZONTAL TAIL. FLIGHT CONDITIONS; AIRSPEED, 135 M/S, A AO A; 0 DEGREES, DROP SIZE; 20 
lim, STATIC TEMPERATURE; -9.3 C, STATIC PRESSURE; 89867 Pa, ICING TIME; 30 MINUTES. LWC; 
.695 G/M 3 . 
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